J. Stat. Phys. 95, May 1999 
ANALYTICITY IN HUBBARD MODELS 

DANIEL UELTSCHI 
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Ecole Polytechnique Federale de Lausann^ 

Abstract. The Hubbard model describes a lattice system of quantum particles with 
local (on-site) interactions. Its free energy is analytic when f3t is small, or fit^ /U is 
small; here, (3 is the inverse temperature, U the on-site repulsion and t the hopping 
coefficient. 

For more general models with Hamiltonian H = V + T where V involves local terms 
only, the free energy is analytic when /3||T|| is small, irrespectively of V . There exists a 
unique Gibbs state showing exponential decay of spatial correlations. These properties 
are rigorously established in this paper. 

Keywords: Hubbard model, local interactions, analyticity of free energy, uniqueness of Gibbs 
states. 



1. Introduction 

Electrons in condensed matter feel an external periodic potential due to the presence of 
atoms. A natural basis for the Hilbert space describing the states of the electrons consists 
in Wannier states, that are indexed by the sites of the lattice. The Hamiltonian for this 
Statistical Physics system can be written in second quantization in terms of Wannier 
states, and with some simplifications we obtain a lattice model [HuH]. 



Forgetting the initial physical motivation, we can consider a lattice model as if the 
particles were really moving on a lattice, and develop a physical intuition in this case. It 
may help in understanding the behaviour of the system. 

The most famous lattice model for the description of quantum particles is the Hubbard 
model. It consists in a hopping term (discretized Laplacian) representing the kinetic 
energy, and a Coulomb interaction between the particles. This interaction is local, or 
on-site, meaning that it is expressed in terms of creation and annihilation operators of 
a same site. Many interesting properties of the Hubbard model have been rigorously 



established, see [ Lieb ] for a review; however, the basic questions about magnetism and 
superconductivity are still unsolved. The present paper brings a modest contribution 
in the sense that interesting phenomena are definitively excluded for some values of the 
thermodynamic parameters. 

The general setting is as follows. We consider a class of models that include the Hubbard 
one, with Hamiltonian 

The vector ^ represents a finite number of parameters, such as chemical potential, mag- 
netic field, ... = 'Ylx&'' ^ of local operators, and T'^ = X^AcZ" ^ 



^Present address: Dept of Mathematics, Ru tgers University, 110 Frelinghuysen Road, Piscataway, New 
Jersey 08854-8019; ueltschi@math.rutgers.edu; http://math.rutgers.edu/~ueltsch] 
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finite-range or exponentially decaying quantum "interaction". The free energy is shown 
to be analytic in /x and (3 in the domain 



r^||e'=l^l < const (1.1) 



A 

where c is a constant depending on the lattice and on the dimension of the local Hilbert 
space. What is quite surprising is that the domain does not depend on the local interaction 
V. The reason is that when T is small with respect to (3, the sites of the lattice are almost 
independent, and the (mean) free energy is essentially that of a model with only one site. 
Such a zero-dimensional system is free from phase transitions, hence its free energy is 
analytic. High temperature expansions would yield comparable results; however they do 
not only require the condition ( |1 . 1\ ) , but also f5\\V\\ < const. 

Concerning the Hubbard model, it is not only true that we have analyticity however 
strong is the repulsive potential; the latter favours this phase, i.e. the stronger the inter- 
actions, the larger the domain. This last result holds at half filling, and if the ratio t/U is 
small enough. In fact, at half filling the Hubbard model is unitarily equivalent to a series 
m powers of t/U that starts with V and the Heisenberg model pSQ| , [MGY| ]: 



-^Hubbard — V + i^Heisenberg + 0{jj^ 



(a rigorous statement can be found in [ DFF|] ). Such a model enters our class, with ||T|| ~ 
jj, hence the condition pt'^/U < const. 

Another example is the Falicov-Kimball model |GM|; it is a Hubbard model where 



only particles of a given spin have hopping, the others being considered as heavy, static 
classical particles. The statements for the Hubbard model are also valid in this case, and 



were proven by Kennedy and Lieb | KL |. 



Section 2 contains precise definitions, statements and proofs for general systems with 
on-site interactions. Section 3 is devoted to the Hubbard model; domains of parameters 
where analyticity can be rigorously proven are proposed with explicit bounds, in the case 
of the 3D square lattice. Finally, the paper ends with a discussion of the Bose-Hubbard 
model, for which partial results may be obtained. 



2. Models with local interactions 

2.1. General framework. Let us be more precise and introduce the mathematical frame- 
work. Let L a i/-dimensional lattice; for instance, L = Z*^, but any other periodic lat- 
tice can be considered. We denote with A a finite subset of L, and the thermodynamic 
limit limA/L/A means lim„^oo /a„ with any sequence of finite volumes (A„) such that 
An £ A„+i, and lim„^oo |c^A„|/|A„| = 0, where 5A„ is the boundary of A„. Let il. a finite 
set with \Q\ = S; we consider the set of "classical configurations" fi^. The Hilbert space 
TC\ at finite volume A is spanned by the classical configurations, i.e. each vector of Ti.\ is 
a linear combination of vectors \nx), n\ G 0^. 

A quantum interaction T is a collection {Ta)a<zi., where Ta is a self-adjoint operator 
with support A. Its action is defined on each T-i\ with A D A, and we have factorization 
properties 

{nx\ TaTa' In^) = (haI Ta |n'^)(nA'| ^a' |n^,) (2.1) 

when j4 n j4' = (a hopping matrix is an example of a "quantum interaction"). Let 
us introduce the connected cardinality \\A\\ of A C L as the cardinality of the smallest 
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connected set containing A, i.e. 

\\A\\ = min \B\; (2.2) 

BdA, connected 

notice that II^H = |^| when A is connected, and ||^|| > \A\ when it is not. We define the 
norm of an interaction to be 

||r||e = supj^||r^||e'=ll^ll (2.3) 

where c is a positive number. Here, ||T^|| is the operator norm of T4. We call V a local 
interaction (or on-site interaction) if I/4 = for all | A| ^ 2; local interactions are denoted 
by (Vx) instead of (V^^y). 

Let /I G M'^ be thermodynamic parameters. The finite volume Hamiltonian Hj^ depends 
on fj, and is given by 

< = E ^-'^ + E ^A- (2-4) 

xeA AcA 

We suppose here that both (Vi^) and (T^) are translation invariant, although periodic 
interactions could be considered with only small modifications. The free energy is given 
by the limit (whenever it exists) 

/(/3,;.) = -ilimi^logTl.e-^''r (2.5) 

We write /o for the "classical free energy" 

/o(/3,M) = -^log ^ {n,\ e-^^" \n,). (2.6) 



We notice that fo is also given by (2^) with T'^ = 0. 

A Gibbs state is a functional that attributes to any bounded local operator K the value 

Tr K e~^^^ 

{K) = lim ^ \„ . (2.7) 

^ ' A/L Tr e-PUt, V ; 

A Gibbs state is exponentially clustering if for any two local operators K and K' there 
exists Ck,k' < oo (with Ck,k' = Ct^K,tyK' for any translations tx and ty, x,y E L) such 
that 

\{KK') - {K){K')\ ^ CK,K'e-'^^^'^">^^ (2.8) 

for some finite constant ^. Here, d{K, K') is the distance between supports of K and K' . 

2.2. Uniqueness of the Gibbs state. The Hamiltonians we consider possess many 
symmetries. For instance, they have translation invariance by assumption; and typical 
models have further conserved quantities, such as the total number of particles, or total 
spin... 

Usually Gibbs states obtained with free boundary conditions ( |2.7D , or with periodic 
ones, have same symmetry properties than Hamiltonians. To obtain pure states with 
symmetry breaking, there are mainly two ways: to introduce boundary conditions, or to 
perturb the system. 

In the quantum case, boundary conditions may be defined by means of a suitable bound- 
ary interaction = (9^)acLi where operators d\ are non zero only for subsets A that 
touch the boundary of A. The corresponding Gibbs state is defined by the expression 
(|23), with Hp, replaced by //a + EacA ^a- 
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Ferromagnetic states are associated with operators of the form (n^.^ — nj-j) apphed on the 
boundary of the volume, while for antiferromagnetism we would use ( — — rixi). Of 
special importance are boundary conditions that break conservation of the total number 
of particles. A state displaying superfluid behaviour should be sensitive to the operator 
^( e~'^ ci+ e'^ Cx) where the sum is over sites touching the boundary. The order parameter 
for superfluidity is Cq (creation operator of a particle at site 0) [P0|, and with above 
boundary conditions we may have {ct)^ = ae'^ with a > 0, revealing the presence of 
superfluidity. 

For a superconductor with Cooper pairs described by a Hubbard-like model, relevant 
boundary conditions are of the form ^( e~'^ c^.|C^| + e^^ CxiCyi) with the sum taken on 
sites of the boundary, close to each other. This allows expectation values of the form 
(cLcL) to be non zero, and this should be the indication of superconductivity |Yang ] . 



The second way to obtain states with less symmetry than the Hamiltonian is to add a 
perturbation, that is then set to zero. As for superfluidity, a good perturbation to consider 
is /i e~'^ cj)+ e^^ Ca;) (see e.g. [ Hua| ) , and the question is whether lim/i^o(c|))^''^ differs 



from zero. 

Here we shall speak of uniqueness of the Gibbs state if it is insensitive to both boundary 
conditions and to external perturbations. In the range of parameters we consider, systems 
are described by Gibbs states sharing the two properties 

• lim\yi^{K)^ does not depend on the boundary conditions d^, provided is 
small enough (independently of A); 

• for all quantum perturbation P with exponential decay, \\P\\c < oo for a large enough 
c, and all local observable K, 

Tr K e"^(^A+a Eaca Pa) 
{K) = lim lim ^ ^ (2.9) 

Notice that P is not necessarily translation invariant, it may even not be periodic. 

Remark: the stability against perturbations can be given a simpler, however more 
abstract definition. Let us consider Q, the Banach space of interactions with finite norm 
(|2.3|), and Q the space of Gibbs states obtained with periodic boundary conditions; G 
is a topological space with the weak topology. Let g denote the corresponding mapping 
Q ^ Q. It is continuous at H £ Q provided g~^{G) is a neighbourhood of if G is a 
neighbourhood of g{H). 

Then the stability of a Gibbs state with respect to perturbations amounts to saying 
that g is continuous at H. 

Indeed, we can see ab absurdo that (|2.9|) implies the continuity of g: suppose G is a 
neighbourhood of g{H) such that g~^{G) is not a neighbourhood of H; since Q is a metric 
space, there exists a sequence (Hn), — > with Hn ^ g~^{G); by (|2.9|) , g{Hn) — > g{H), 
then g{Hn) € G for n sufficiently large, and therefore G g^^{G). Conversely, for any 
open set G that contains g{H), g~^{G) is a neighbourhood of then if H, we 

have Hn G g^^{G) for n sufficiently large, therefore g{Hn) G G. 

2.3. Result. In order to state the result, we let H be a constant that depends only on the 
lattice, such that 

#{A 9 X, connected, \A\ = /c) ^ 2^. (2.10) 
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A possible choice, probably not optimal, is 13 = (2z^)^ for the z^-dimensional square lattice. 
The Golden Ratio appears here, that we write cj) = following a standard convention.0 

Theorem 2.1 (Analyticity in models with local interactions). 

Assume that and T^^ are smooth, i.e. that all matrix elements of Vi^ and are 
analytic in fj, for all x and all A. Let c ^ cq = log S + log 23 + (p + 2 log (p. Then in the 
domain 

/^r'^iic < 1, 

(i) the free energy exists in the thermodynamic limit and is analytic in (3 and fx; 

(ii) the Gibbs state converges weakly in the thermodynamic limit; 

(iii) the Gibbs state is exponentially clustering with a correlation length bounded by = 
4(c-co)-i. 

And the Gibbs state is unique, i.e. 

(iv) the Gibbs state is stable with respect to boundary conditions with /?||T'^ + 9^||c < 1 
for all A; 

(v) the Gibbs state is stable with respect to all external perturbations P with \\P\\c < oo. 

Remark: the bound 4(c — cq)^^ for the correlation length is rather arbitrary and could 
certainly be improved. 

The stability against boundary conditions should hold for any bounded d^, not only 
small ones. However, such a statement is difficult to prove in quantum systems, where we 
have to deal with negative or complex numbers. 

Proof of Theorem \2. j| (i). The idea of the proof is to expand the operator e~^^^ with 
Duhamel formula; it allows next to express the partition function as the one of a polymer 
model. After having shown that the weights of polymers have exponential decay with 
respect to their size, the analyticity of the free energy is a result of cluster expansions. 
The Duhamel formula (very similar to the Trotter formula) yields 

Tr e-^^A = TV g-^^-eAKT + ^ ^ f ^ 



m^l Ai,...,A„iCA 



0<Tl<...<Tm<f3 



IV e'^^^^eA^^" e'^^^^^^^^^eA^^" . . .T^ Q-W-'rni)'ExeA^^ _ (2.11) 

For given Ai, . . . , Am, we construct the graph Q of m vertices, with an edge between 
i and j whenever Ai f] Aj 7^ 0. Decomposing Q into connected subgraphs, it induces a 
partition of {Ai, . . . , Am} into i subsets {i ^ m). We let ^1, . . . , C to be the 
unions of sets Ai, . . . , Am for each partition. As a result, to each sequence Ai, . . . , Am 
corresponds a unique set {^1, . . . ,Ai} of subsets of Z,'^, such that 



(2.12) 



We call Ai, . . . ,Ae polymers and define their weight 

p{A) = e^M^'^^^\^\ Y: i-ir E E / dn...drm 

(n^l e-^i^-e^^^^^T^ e-^r2-r,)Y:.eAy-^ . . . g-^^"""™) ^-e^ |n^). (2.13) 



^It is a pleasure to welcome here the Golden Ratio. Its presence is however fortuitous and does not 
involve any of its special and beautiful properties. 
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The sum is over sets Ai, . . . , Am satisfying two restrictions: (i) U'^-^^Ai = A, (ii) the graph 
Q defined above is connected. The partition function can then be rewritten as 

TV e-^^A = e-^/o^/^'/^)!^! HpCA). (2.14) 

Mi,...,A}i=i 

We have now to bound p{A); first the matrix element: 

m 

i=i 

Let Cq be the lowest eigenvalue of Vi^; since /o(/3, /i) ^ e^, we have 

||e-/3^"|| = e-^^o ^ ^-PMP,f^)_ (2.16) 

Furthermore \^-^\ = S*'"^' and the integral over "times" {tj} brings a factor l3"^/m\; using 
^ Z^JLi ll^jll) obtain 

fom ™' 

\piA)\ ^ 5l^le-ll^ll E ^ E nil^A.II^^"^^" 

m ^ 1 CAj=l 

^ g-(c-iogS-i)Mll _ (2.17) 



Results of cluster expansions are summarized in Proposition |2.2| below. From this we 
obtain the following expression for the free energy 

It does not depend on x, because the Hamiltonian is translation invariant. Since $"'"(C) 
is analytic in P and /x, and the series converges uniformly, the free energy f{P,fj,) is an 
analytic function by Vitali theorem. □ 

Proposition 2.2 (Cluster expansions). 

Let us recall that we can choose U = (2i/)^ for the u- dimensional square lattice and 
4> = ^2'^ is the Golden Ratio. 

Assume that a function z^'^ : 'P(L) — > C is given and such that for all Ad'h, 

• \zf^'t^{A)\ ^ e-^ll^ll withr ^ tq = log 2n + - 1 + 2 log 

• z^'^{A) is analytic in l3,fi. 

Then there exists an analytic function : V{V{1L)) —>■ C such that 

log Y: fl-'''^(^.)= E '^^(C); 
{Ai,...,Afe} j=l C={Ai,...,Ak} 
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Let supp C = Uasc^;' ^"^{C) = if C is not a cluster, i.e. if C = CiU C2 with supp Ci H 
suppC2 = 0. This function has exponential decay: 



C={Au-,Ak} 



1 



for all X Gh. Here, we set \\C\\ = X^^gc* ll^ll- 

This proposition is an immediate corollary of Kotecky and Preiss theorem on cluster 



expansions |KP(| (see |Dob| for an elegant and simpler proof). To make the link between 
our notation and theirs: 



Here 




A 


7 


23 


K 




07 


(r-ro)PII 


d(7) 



Our polymers are not necessarily connected; but we note that their entropy satisfies 
#{A B X, \\A\\ = A;) ^ 2'' #{A 3 x, connected, \A\ = k) 

Hence the factor 2 in front of 

A useful consequence of this proposition is the existence of the thermodynamic limit of 
the free energy of a gas of polymers; if the weight z^'^ is periodic with respect to lattice 
translations, then the limit 

^ ' ' {A^,...,Ak} J=l 

yliCA,AinAj=0 

exists and is analytic in (3 and /x. 



Proof of Theorem \2. j| (ii). Having the expansion ( 2.14| ) for the partition function, the 
treatment of expectation values of local operators is standard. The expectation value 
of K involves the quantity Tr K e~^^^ , that we expand as before with Duhamel formula. 



TrKe-'^^A =TVA'e-^S^sA^" + ^ (-1)" E / 



dri 



, dr„ 



yli,...,A„CA-^°<^i<-<^™<^ 



e 



(2.19) 



We construct the graph ^ of (m + 1) vertices, for the sets suppET, Ai, . . . , Am, and we 
look for (l + 1) connected components. One of these components contains suppii', and 
we denote it by Ak', others are denoted by ^1, . . . ,Ai, as before. The weight of Ak is 
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modified by the operator K; namely, 



Pj,[Ak) = e^^«(^''^)l-^^^l jTri^e-^^^^^K^" + ^ (-1)- ^ 

j dn . . . dT„(n^^ I e"^^ ^-^^K T^^ e"^"^^""^^ ^^^^k . . 

J0<Tl<...<Tm<l3 ^ 

___T^^^-(P-r,n)i:.,^,v.^ I^^k)}- (2.20) 

The expectation value of K takes form 

<^'>=i7Lzi(b) s «(-^'^-)nf(A) (2.21) 

where the sum is over non intersecting sets Ak,-Ai, . . . ,A£, i 0. The weight pK has 
also exponential decay: 

\pk{Ak)\ ^ ||i^||e^!l«"PP-^"ll e-("-'°s-^-^)"-^^^''l . (2.22) 

From cluster expansions, we obtain an expression for {K): 

W=|im Yl PK{AK)exp{- Yl ^^(C-)}. (2.23) 

AkCA CsuppCcA 
Ak^suppK suppCnAK¥=9 

The limit exists, because the sums converge uniformly in the volume A. This proves the 
weak convergence of the Gibbs state in the thermodynamic limit. 

□ 



Proof of Theorem \2.l\ (in). Using cluster expansion, this is standard stuff. Expanding the 
expectation value {KK') as before, we get 

{KK')= Yl PKK'{AKK')exp{- Yl ^^(^)} 

Akk' 3supp i^Usupp K' C,supp CnAjfi^i ^0 

+ Y PK{AK)PK'{AK')exp[- Y ^''(C)}- (2-24) 

AkDsuppK C,snppCn[AK^AK,]^<l> 
AifiDsupp K' 

AKnAi^/=iIi 



The weights pK , Pk' , Pkk' are given by ( 2.2C| ) when substituting K with K, K' , KK' 
respectively. 

We define {KK'f'^'''^ to be as (p|), but with sums only over polymers and clusters 
of connected cardinality less than ^d{K,K') (we call such polymers short, they are big 
otherwise). We denote {KK')^'^ = {KK') - {KK'y^°'K 

When the expectation values are restricted to short polymers and clusters, correlation 
functions are zero: 

{KK'Y^"'^ = (if)'^hort ^^/^short ^2.25) 

Therefore 

{KK') - {K){K') = {KK')^'^ - (i^)^'s {K') - (K)«hort ^_^/^big ^ ^2.26) 
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Expectation values involve sums over polymers of connected cardinality bigger than 

^d{K, K'), and this has exponential decay. There are also terms 

exp{- y: ^^(c^)}-exp{- y: ^^(C)} 

|suppC| s; \d{K,K') 

= exp{- J2 ^^(C)}[exp{- ^''{C)}-l\. (2.27) 

C,supp cnAj^d C.supp cnA^d 

|suppC| ^ \d(K,K') \svLppC\>\d{K,K') 



We know from Proposition 2.2 that the sum over clusters has exponential decay; more 
precisely, the quantity between brackets is bounded by C{A) e~*^^^'^ )(c-co) ^ -y^here C{A) 
depends on |^| only. 

Exponential clustering is now clear. 

□ 



Proof of Theorem \2. l\ (iv). The proof is rather standard, so we content ourselves by out- 
lining it. 

Expanding Tr e~I^^A~/^^AcA9A with Duhamel formula, we obtain an expression very 
similar to ( |2.19| ). The difference is that operators now appear in the second line of 
(|2T9| ). We define pk{Ak) to be as in (l2.2C| ), except that at least one boundary operator 
shows up. Similarly, let p{A) be like ( |2.13 ) but with at least one d^; an important 
property of these weights is that they are zero if the polymer does not touch the boundary 
of A. As a result we get 

= E {pKiAK)+PKiAK))exp{- Y {^''iC) + ^^iC))} (2.28) 

AkCA CsuppCcA 
Ax^suppK suppCnAK¥=9 

where ^"^{C) is such that [<I>'^(C) + <I'^(C)] is the truncated function for polymers with 
weights [pIa) + p{A)]. 

A few more developments lead to an expression for (K)^ that is equal to {K)\, plus 
terms that connect supp-ftT with the boundary of A, and that decay exponentially quickly. 
In the thermodynamic limit, this correction vanishes, and therefore for all d^: 

|im (i^)f = {K). (2.29) 

□ 



Proof of Theorem 2A (v). Once we have an expansion in terms of clusters, the proof of 
stability with respect to external perturbations is actually fairly simple. 

First we note that if Ylm ^o\9n\ < oo uniformly in a, and g-a when a ^ 0, then 

Y^n^Y 9n- 

n ^ n ^ 

Indeed, for any e > there exists m such that 

e 



E {\9n\ + \9n\) ^ 



n > m 
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Furthermore, for all n there exists a„ > such that \g!^ — 9n\ ^ 2m '^^^n a ^ q„. 
Choosing a = min„<m ctn, we have for all a ^ a 

\j2 9n-Y.9n\ ^ ^\9n- 9n\ + Yl (l^nl + M ^ E. 
n ^ n ^ n<m n ^ m 

We add to the Hamiltonian a new quantum interaction P with ||P||c < oo. Since 
||T'^||c < 1, there exists q > such that for all a ^ a, 

||r'^ + aP||c < 1. 
Retracing the steps above, we obtain 

(K)"= Yl pU^K)exp{- Yl '^"JC)^ (2.31) 



(2.30) 



where pf^{AK) is obtained by replacing with (T^ + oPa) in (|2.2C|); s imilarly, <I>^ is 
constructed with weights that we get by the same substitution in ( ^.13 ). 

The weights and p'^ are absolutely convergent series in matrix elements of {T^ + 
oPa}, therefore p"' ^ p and p'^ ^ px as a ^ 0. Hence —>■ and also 



c,suppCn^A'^0 c,suppCnylK^0 

Since the expression ( |2.31| ) for (i^)" is absolutely convergent, this implies that {K 
{K). 



2.32) 



□ 



3. The Hubbard model 

The phase space of the Hubbard model is the Fock space of antisymmetric wave func- 
tions on A C L. A convenient basis is the one in occupation numbers of position operators. 
It is given by { |nA)}„^gQA where 0, = {0,t,|,2}. The Hamiltonian is 



HA = -t Y 4<TC?;a + ^^E"^-t"^i~/^E("^t + '^xi) 



(3.1) 



<c,y>CA 

^e{T,i} 



xeA 



zeA 



where the first sum is over nearest neighbours x,y € A. The first term represents the 
kinetic energy, the second one is the local repulsion between electrons, and the last term 
is the chemical potential multiplying the total number of electrons. 



(a) (b) (c) 

Figure 1. Domains of analyticity, stemming from (a) high temperature 



expansions, (b) domain T>i of Theorem 3.1, and (c) domain T>2 of the same 
theorem. 

High temperature expansions yield analyticity of the free energy for all /? such that 

Pt < const and PU < const. 
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see Fig. 1; the domain of analyticity may be extended, as we see now. Let x be the 
maximum coordination number of L (x = 2i/ for z^-dimensional square lattice), and e 



a2^-1 



Theorem 3.1 (Analyticity in the Hubbard model). 

Let A = min(//, U — Thermodynamic limits of the free energy and of the Gibbs state 
exist in the domain V1UV2, where 



Vi = {{(3,t,fi) ■.(3t<e} 



and 



V2 = {{P,t,n) :0 < fi<U and < 2xe^il - 2t/eA)}. 
The free energy is analytic, and the Gibbs state is exponentially clustering and unique 



(that is, stable against boundary conditions , 
\\P\\c < 00 for a sufficiently big c). 



< 1, and external perturbations P, 



1 



Remark that the domain P2 is meaningful only if t is small enough, namely ^ < ^e. 
For the 3-dimensional square lattice, we find for T>i the condition [it < 1.75... • 10~^, and 
for 1^2, ^ < 3.68... • 10~^(1 — 1.14... • 10^-^). Of course, the domain of analyticity is much 
larger than these domains, where analyticity is proven to hold. 

These properties likely hold for all /3 < 00 in dimension 1, and possibly also in dimension 
2. When v ^ 3, a domain with antiferromagnetic phase is expected for t <^ U and 
(5t^ /U > const. Such a phase can be proven in the asymmetric Hubbard model, where 
electrons of different spins are assumed to have different hopping parameters [ [KL| , LM, 
pVIM] , pFF| ] (see also [DFFR] and [KU] for two general methods to study rigorously such 
situations). Assuming this to be true in the standard Hubbard model, we observe that 
the condition for the domain P2 is qualitatively correct; this is illustrated in Fig. 0. 




temperature 



uniqueness 



antiferromagnetic 



Figure 2 . Phase diagram of the Hubbard model. Antiferromagnetic phase 
is expected for dimension z/ ^ 3; it can be proven when z/ ^ 2 for the 
asymmetric model. 



Proof of Theorem 3J_, domain T>i. The classical free energy of the Hubbard model is eas- 
ily computed and is given by 

/o(/5,/i) = -;^log[l + 2e^^ + e-^^+2/3;.]_ (32) 

Let us write the kinetic operator —t Ta where A = {<x,y>, a) and the notation 

A C A means x, y G A; Ta = cLaCya- 
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The expression (|2.13| ) for p{A) takes the following form 



p{A) 



m 1 



dri . . . dr„ 



Ai,...,Am n^eQ-^ •^°<^i< 



.<Tm</3 



with a restriction on the sum over Ai , . . . , , namely their union yields A and they are 
connected in the sense of the graph G described above. Notice that here p{A) = if A is 
not connected. The expression for Pk{Ak) is similar, compare with ( ^30|) . 

We obtain the domain Di of Theorem |3.l| by proceeding as before. Namely, we bound 
the matrix element with e"^*^"!-^' ^ Q-Pfo{P^i^)\-A\ ^ observations allow to slightly 

optimize the bound for p{A). First, there are no more than '^\A\ sets of nearest neighbours 
in A. Second, if we choose Ai, . . . , (such that they cover A), then there is at most 
one configuration such that (n_4| Ta^ ■ ■ ■ Ta^ \n_A) differs from 0. 

Therefore we obtain the bound 

ipty 



\p{A)\ ^ e 



-{c+l)\A\ 



^ e 



-c\A\ 



assuming that 



E 

m > 



2xe^+i/3t ^ 1. 



-(— |A|)™'4™' e*-'^^"^'*'" 



(3.4) 



(3.5) 



Here the polymers are connected sets; in this case. Proposition 2.2 holds with 2^ replaced 

by n. 



Therefore the condition ( |3.5| ) must be fulfilled with cq = logD + — 1 + 21og0. 
When this inequality is strict, it also holds with c > cq, so that we obtain exponential 
clustering and stability against perturbations or boundary interactions. 

□ 

Proof of Theorem \3. 4 domain T)2- Domain P2 benefits from the following geometric rep- 



resentation (see Fig. for intuition). First we let = n^, then 



m— 1\ 
A 



±Ta 



\n 



A'^ 



\n 



A/ 



±TAn \n 



A/^ 



\n 



A/ 



^Tai \n^). The last condition follows by cyclicity of the 



trace. 




Figure 3. Three loops. 
Next, if Aj = {< Xj,yj >, cr), we define horizontal bonds B C M.'^ x [0, (3] per 

m 

B = U XjUj X {r, } 





1 T 2 
2 

1 T 

2 

2 I 
I 



2 i 
2 

i 

T 



2 T i -I 2 t 



A 
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where xji/j C M'^ is the segment joining xj and yj. We consider vertical segments 

S = U^{x G L : ni G {0, 2}} x [r,, r,+i], 

where we set tq = 0, r^+i = (3 and nP = n™. Actually, 5 is a subset of L x [0, /3] per; but 
with a small abuse of notation, we consider 5 C M*^ x [0, /3]per- 

The set BuS decomposes into a finite number of closed circuits that we call /oops.^ To 
be precise, a loop £ is a pair (supp£, where supp£ C M"^ x [0, /3]per is the support 

of £, and A{£) = (Ai, . . . , A^^^)) are successive applications of operators Taj , • • • , Ta^^^^ 
here m{£) is the number of horizontal segments ("jumps") in £ that we mark out 1, . . . ,m 
in increasing vertical coordinates, and Ta is the operator associated with the segment j. 

The weight p{A) can be written as an integral over sets of loops, with many restrictions. 
In particular, each vertical line {x} x [0, /3]per, x £ A, must intersect at least one loop. As 
a consequence, to a given set of loops corresponds at most one sequence of configurations 

p(^) = e'^^''^^'^)!^! e'^^l-^l ^ / d^i . . . d4 e(^i, • • • ,4) H ^(^i) (3-6) 



k > 1 



where 



and 



Ae(£i,...,4) 

z[^e)=f^{e)^-M.-{v-m^ ^ (3.8) 

Let us explain these notations. The configuration n_4 is defined by (^i, . . . ^£k)\ namely, 
if {x} X {0} G supp^j, we know that Ux G {0,2}; looking at the first occurrence of an 
operator T^, A 9 x, we can check whether a particle is created or annihilated at x, in 
which case = or = 2 respectively. Similarly, if {x} x {0} ^ Ujsupp£j, we have 
"^x £ {T) i}; if the first operator Ta such that A 9 x creates an | electron, or annihilate a 
1 electron, we have nx=i\ otherwise nx=\. 

The product is over all operators Ta that occur in the loops, ordered in decreasing 
vertical coordinate of the corresponding horizontal segment. Notice that e(-) G {—1, 0, 1}. 

The vertical length of a loop is \£\ = \£\o + \£\2, where \£\j denotes the length of all 
vertical segments where the configuration takes value j. We bound 

z{£) ^ t'-W e-^l^l . (3.9) 

We have the following bound for p{A) (we use /o(/3,/u) ^ /x): 

\p{A)\ ^ e-(^+i)l-^l d^te^+M'^^'^e-^l^lj". (3.10) 



k > 1 



The integral over one loop with m jumps may be evaluated in the following way. 

1. We choose two nearest neighbour sites in A (there are less than x\-^\ possibilities), 
we integrate over a number r in [0, and we choose a spin; we obtain the first jump 
of the loop. 



^This representation has many similarities with that of [MM|, introduced for the Fahcov-Kimball modeL 
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2. We decide whether the loop is going up or down in the vertical dimension, we in- 
tegrate over the vertical distance, we choose a neighbour of our site and a spin; 
integration over the vertical distance is bounded by 

1 

A- 

We repeat this procedure until the last jump but one. 

3. For the last jump, we decide whether the vertical direction is up or down, and we 
integrate over the distance, yielding a factor 1/A. Then the loop completes itself 
in a unique way (provided there is a way); therefore there are no sums over nearest 
neighbour and spin. 

Since the first jump is arbitrary, we can divide by m the contribution of loops with m 
jumps. Notice that the second step is superfiuous when m = 2. We obtain 

Y [ ,,,,,, d^(te'=+i)™e-^l^l ^ y -{te^+^rx\A(32[2^x^r~^2^ 



m^2 m{i)=m m ^ 2 



^ |^|^2xe2(^+^)— -i^. (3.11) 
A 1 - 4x e^+i ^ 

From the conditions s < 4^ e~('^+^) and ^ < ^ e~^(^+^) (1 - 4x6^^+^ ■^), we finally have 

\piA)\ ^ e-(^+i)l-^l rfl^l'' < e-^l-^l (3.12) 

with c = log!3 + (p — 1 + 21og(/>. Since the weights of polymers are analytic functions of 
/3, 11, so is the free energy in the thermodynamic limit. 

Existence and properties of the Gibbs state are readily obtained by repeating the proofs 
of Section 2, using above estimates. 

□ 

4. The Bose-Hubbard model 

The Bose-Hubbard model describes a lattice system of interacting bosons. In a finite 
volume A G L, the phase space is the Hilbert space with basis { I^^a) : riA G n^}; the 
Hamiltonian consists in a kinetic operator and a local repulsive interaction: 

HA = t J2 4cy + uY,if^l-n^)- f^^n^- (4.1) 

<x,y>cA- xeA xeA 

The first term is a standard hopping operator between nearest-neighbours; the second 
term describes the local repulsion between bosons (each pair of particles at a given site 
contributes for 2U); the chemical potential /i G M controls the density of the system. 
It has been introduced in pWGF | and despite its simplicity, it has very interesting 



phase diagram, see Fig. ^ A phase transition insulator-superfluid is expected when the 
hopping coefficient increases. Of course one would like to have mathematical statements 
to support this, but the superfluid phase of interacting particles is hard to study.^ On the 
other hand, one can tame the insulating phase much more easily. When t/U is small, and 
2U{k — 1) < /X < 2Uk, it is possible to show that the Gibbs state exists at low temperature 



4 



The best rigorous statements concern the hard-core Bose-Hubbard model, where off-diagonal long- 



range order can be proven using reflection positivity for special value of /i [ DLS ; this constitutes a beautiful 
result, although it does not allow to study pure states. 
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and 4 e is close to the projector onto the configuration Ux = k for ah x G L; moreover, 



the density of the (quantum) ground state is not only close, but equal to k |BKU 



= 4 



6U 



p = 3 



4(7 



2U 



Figure 4. Zero temperature phase diagram for the Bose-Hubbard model. 
Lobes are incompressible phases with integer densities. 



superfluid 



We show here that these phases can be reached from the high temperatures without 



phase transition (Fig. |5|). Theorem 2.1 does not apply here, because the single site phase 
space is infinite. Actually, boson systems and unbounded spin systems present some 
difficulties at high temperatures, since partition functions diverge at /3 = 0. Results for 



small P have been obtained in |PY|. We can prove analyticity of the free energy and 



existence of Gibbs state when (3t is small, but we are unable to show stability against 
perturbations, or against boundary interactions, that do not conserve the total number of 
particles (see discussion in Section ^ ). 

Because the phase space has infinite dimension, we need to define the Gibbs state as a 
functional over possibly unbounded operators, as for instance ct, or number operators; but 
not all operators can be considered. In order to define a suitable class of local operators, 
let Na_ be the number operator in A with minimum eigenvalue 1, i.e. 



Na \nA) 



pa; 



if = for all x £ A 



\nA) otherwise; 



(4.2) 



this operator has an inverse which is defined everywhere. Defining the boson norm of a 
local operator K by 



\K 



I boson 



sup 



{n\ iV^ 



supp.ft'-^^suppis: \^ 



(4.3) 



we consider the class /C of local operators with finite boson norm. It is not hard to check 
that 1 1 TV 



^ 1, and thus CxCy G /C. 



Theorem 4.1 (Analyticity in the Bose-Hubbard model). 

There exists a function e{U, /i) > such that for all j3 with f3t < e, 

• the free energy f{P,fj,) exists in the thermodynamic limit and is analytic in (3 and fi; 

• the Gibbs state {■) : fC ^ C, with free or periodic boundary conditions, converges 
weakly in the thermodynamic limit and is exponentially clustering. 

Remark: adding a hard-core condition on the model, e.g. by limiting the number of 



bosons at a given site to M < oo, Theorem 4.1 becomes a consequence of Theorem 2.1 
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(with moreover uniqueness, and absence of superfluidity). Actually, we expect Theorem 



to hold when T has finite boson interaction norm 

jbosons ^ sup^ ||iV^^TAiV~^||e^il^ll 



(4.4) 



However, we are unable to handle such a general situation; the Bose-Hubbard model 
conserves the total number of particles, and this property plays a crucial technical role. 



The present result, together with [ BKU ], shows that the free energy is analytic in a 
domain that includes low and high temperatures, corresponding to insulating phase; see 
Fig. |. 




1 





uniqueness 



Figure 5 . Analyticity of the free energy and existence of the Gibbs state 
can be proven on the left of the grey frontier. 



Proof of Theorem |^.4 Expanding Tr e~^^'^ , we obtain (2.20) and ( p. 21 ). We show now 
exponential decay for pxi'A.K)] similar (and simpler) considerations lead to exponential 
decay for p{A). Analyticity of the free energy and existence of Gibbs state are then 



consequences of Proposition 2.2 for cluster expansion. 
The condition Uj^iAj U suppii' = Ar implies 

\Ak\ ^ \suppK\+m 
(recall that \Aj\ = 2 for all 1 ^ j ^ m). Then 



(4.5) 



\pk{Ak)\ ^ e-^l-^^l e^l""PP^I e'^Ml3,f^)\AK\ 



m > 



E 



<Xl,yi>,...,<Xm.,ym>(lAK 



1 

J0<Tx<-<Tm<l3 



dri 



\nAi. 



(4.6) 



For given |n_4^, ), let n\^, ... , nj^ G N-^^ such that 
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Then the matrix element in the above equation takes form 
{nAAK \nX,) e'^ ^^--k^--k\^^ ^^-k^ (n^J 4^c,, \n\^) . . . 

m+l m 

^^eylK i=i j=i 

in the last line we set tq = 0, Tm+i = /? and = nj^j^. We used the inequality 

(n'l c^Cy |n) ^ n^, + n^^. Now 

m+l . . m+l 

i=i j=i 

The sum over n_4^, can be replaced by a sum over f^-;^^; since this last bound does not 
depend any more on , z 7^ j, we can sum over <Xi^yi><Z Ak, and since ^^gAk ~ 
J2x£Ak ^'^^ ^' ^^^^ product of (^3) is bounded by 

m 

n E ^ U E -^)'"- (4-9) 

i=l <ri,S/i>C.4/f xeAx 

Collecting these estimates, we get 

\Pk{Ak)\ ^ q-^\-^k\ ^^J^^^boson ^c\suppK\ ^f3fo{P,fi)\AK\ 

^ (m + l)re'='"^x'"E(l+ E n,)™+^-''^-s^K<"-4xlV'.l-^K). (4.10) 
When (3t is small, we have {f3te'^x)"^ ^ e'^ x)™"*""^! therefore 



E ^^[(1+ E nx)^/:^^]" < exp{V/3ie-+2^(l+ J2 ^-)}- (4-11) 

m ^ xe^x xgAk 

We obtain finally 

\Ak\ 



n=0 



(4.12) 



The sum over n is absolutely convergent, so that the quantity between brackets converge to 
1 when pt — > 0. Having chosen c sufficiently large to ensure validity of cluster expansion, 
we see that pK decays exponentially when pt is small enough. 

□ 
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